Amplitude-Phase Coupling in a Spin- Torque Nano-Oscillator 
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The spin-torque nano-oscillator in the presence of thermal fluctuation is described by the normal 
form of the Hopf bifurcation with an additive white noise. By the application of the reduction 
method, the amplitude-phase coupling factor, which has a significant effect on the power spectrum of 
the spin-torque nano-oscillator, is calculated from the Landau-Lifshitz-Gilbert-Slonczewski equation 
with the nonlinear Gilbert damping. The amplitude-phase coupling factor exhibits a large variation 
depending on in-plane anisotropy under the practical external fields. 



When a direct current / flows into a magnetoresistive 
(MR) device, a stationary magnetic state becomes un- 
stable and a steady magnetic oscillation is excited by the 
spin-transfer torque. The oscillation is expected to be 
appUcable to a nanoscale microwave source, i.e., the spin- 
torque nano-oscillator (STNO).— According to the the- 
ory based on the spin-wave Hamiltonian formalism,— i^i^ 
the frequency nonlinearity plays a key role in determin- 
ing the behavior of the oscillator. It has been shown that 
the strong frequency nonlinearity leads to significant ef- 
fects on the power spectrum of STNO in the presence 
of thermal fluctuation: a linewidth enhancement^ and 
non-Lorentzian lineshapesS. In this paper, the impor- 
tant nonlinearity is examined. From the Landau-Lifshitz- 
Gilbert-Slonczewski (LLGS) equation as the model of 
STNO, we calculate explicitly the magnitude of the quan- 
tity corresponding to the normalized frequency nonlin- 
earity N/Tcfi (see, e.g., Eq. (4) in Ref. [g) of the spin- 
wave approach. In particular, we take account of in-plane 
anisotropy of a magnetic film which has been neglected 
in the early studies^i^i^i^, finding the large effect of the 
anisotropy on the nonlinearity. 

We describe STNO by a generic oscillator model. It 
is known that small-amplitude oscillations near the Hopf 
bifurcation point are generally governed by the simple 
evolution equation for a complex variable W{t) known 
as the Stuart-Landau (SL) equationji The SL equation 
is derived as a normal form of the supercritical Hopf bi- 
furcation from the general system of ordinary differen- 
tial equations. Accordingly, the LLGS equation similarly 
reduces to the SL equation in the case where the Hopf 
bifurcation, which represents a generation of magnetic os- 
cillations in STNO, occurs. The reduction of the LLGS 
equation can be executed by the reduction method based 
on the center-manifold theorem. At finite temperature, 
there exists inevitable thermal magnetization ffuctuation 
in STNO.^'^ We include the thermal effect into the mag- 
netization dynamics by just adding white noise term to 
the SL equation, i.e., STNO in the presence of thermal 
fluctuation is described by the 'noisy' Hopf normal form: 
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tor M (see Eq. (O below). In Eq. (P), f2 represents a 
fundamental frequency, t is a normalized dimensionless 
time, and r](t) is the zero-mean, white Gaussian noise 
with the only non-vanishing second moment given by 
{'ri{t)v{t')) = 4i5(t — i'). p is the bifurcation parame- 
ter. An oscillation is generated when p becomes positive. 
In the context of STNO, p oc [I — Ic) where Ic is the 
threshold current. The parameter 5 quantifies the cou- 
pling between the amplitude and phase fluctuations and 
is called the amplitude- phase coupling factor. It is 5 that 
we calculate numerically in this paper and that corre- 
sponds to the normalized frequency nonlinearity N /TcS 
of the spin-wave approach. The amplitude-phase cou- 
pling factor 5 affects the power spectrum of an oscillator 
and leads to linewidth enhancement and non-Lorentzian 
lineshapesJ^iii Due to its effect, the factor 5 is also called 
the linewidth enhancement factor}^ Eq. ([T]) is often used 
as the simplest model of a noisy auto-oscillator in many 
fields, for example, electrical engineering, chemical reac- 
tions, optics, biology, and so on.^°'l? Therefore, we can 
easily compare STNO with conventional oscillators and 
clarify its features. 

The amplitude-phase coupling factor i5 is obtained in 
the procedure of the reduction of the LLGS equation. In 
the following, we first explain the LLGS equation. Then, 
following Kuramoto's monograph^., we consider an insta- 
bility of a steady solution and execute the reduction of 
the LLGS equation. 

The magnetic energy density of the free layer of STNO 
is assumed to have the form 



M2 



{M -xf + ]-^-kM -N ■ M, (2) 



where Mg is the saturation magnetization, i?oxt = H^x-h 
Hyij -\- HzZ is an external field, is uniaxial anisotropy 
along the x direction, and J\f is the demagnetizing ten- 
sor; N — dia,g{Nx, Ny, Nz). Using the spherical coordi- 
nate system (see Fig. [1} , we describe the magnetization 
dynamics of STNO by the LLGS equation 



cos ?/; = -a{^)ip - Fi (0, ip,ujj) 
-0 = a(0 cos ^ -f F2 {(j), tlJ,ujj) 
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where W is the normalized complex variable represent- 
ing the amplitude and phase of a magnetization vec- 



where Fi{(j),4'i^j) = {'j / Ms)d£ / dip — a{(j),ujj) and 
F2(0,-0,wj) = {-i / {MsCOS^))d£ / dcf) + b{(t),'i}},ujj). 7 is 
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FIG. 1: The spherical coordinate system {(j>,'4') for the direc- 
tion of the free layer magnetization m = M/Ms of STNO. 
p denotes the direction of the pinned layer magnetization; 
p = (cos i/)p cos (^p , cos -i/jp sin (j)p , sin tpp ) . 



the gyromagnetic ratio. The second terms of Fi re- 
sult from the Slonczewski term Tj = {'yaj/ Ms)M x 
(Af X p) in which aj is proportional to the cur- 
rent density J through the free layer^^. Therefore, 
a{4>,ujj) = w J cos ■i/'p sin((/) — and 6(0, -0, wj) = 
ti;j[cos'!/'p sin0cos(0 — 0p) — sin ^p cos i/'] , where loj = 
"fttj. Q;(^)-terms of Eqs. ([3]) are the generalized Gilbert 
damping terms proposed by Tiberkevich and Slavin.— 
We take into account only the first non-trivial term of 
the Taylor series expansion for q;(^) by the magneti- 
zation change rate ^ = (9m/9t)^/(747rMs)^; a(^) = 
Q!g(1 -I- qi£,). According to Ref. il5j, the nonlinear LLGS 
model with qi — Z gives a good agreement with the ex- 
perimental results of Ref. 1 and Ref. T^. 

An instability of a steady solution of Eq. ([3]) is consid- 
ered. A steady solution (0o(wj), '(/'o(wj)) is derived from 
Fi{(j)o, i/jo, wj) = 0. Shifting the variables as ui = — 0o 
and U2 = V' ~ "00 J we have the Taylor series of Eq. ([3]) as 
follows, 



u = Lu + N2UU + N3UUU 
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where u — (1/1,1*2)"'". Here, the diadic and triadic nota- 
tions^ have been used. The stability of a steady solution 
is determined by the eigenvalues of the linear coefhcient 
matrix L: A± = T ± (T^ - detL)!/^. T is defined as 
r = r{ujj) = (l/2)trL and plays the role as a control pa- 
rameter since it depends on uj. We confine ourselves to 
the case where the Hopf bifurcation occurs. Then, A± is a 
pair of complex-conjugate eigenvalues. The point, F = 0, 
is the Hopf bifurcation point; while a steady solution re- 
mains stable for F < 0, it becomes unstable for F > 0. 
The bifurcation point corresponds to the threshold ujj 
which is determined by trL = and i^i(0Oi "^Oi ^j) = 0- 
Near the bifurcation point, we divide L into the two 
parts; L — Lq + FLi, where Lq is the critical part and 
FLi is the remaining part. Corresponding to L, A+ is also 
divided into the two parts; A+ = Aq -l-FAi. Although Li 
and Ai generally depend on F further, we neglect their 
dependence and evaluate them by the values at F = 0. 
Accordingly, Aq — iujQ and 
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FIG. 2: (Color online) (a) Power P divided by RqI^ with Ro = 
13.6 n and (b) linewidth (FWHM) of the signal of STNO as 
a function of applied current /. Dots are experimental data 
at T = 150 K taken from Ref. 'iS. Red lines are theoretical 
fitting curves based on the model of Eq. Q. 



Lq corresponding to the eigenvalue Aq are denoted as U 
and U* , respectively. These are normalized as U*U = 
U*tj — 1 where TJ means a complex conjugate of U . 

Let us apply the reduction method to Eq. (|4]). The SL 
equation for a complex amplitude W{t), 

W ^T\iW - g\W\'^W (6) 
and the neutral solution for the magnetization dynamics. 



(7) 



are obtained within the lowest order approximation.^ Un- 
der the approximation, only the Taylor expansion coef- 
ficients up to the third order are needed. The complex 
constant g in Eq. ^ is given by 



g=vi+W2^-i{U\ N3UUU) 

+ AiU*,N2UVo) + 2{U*,N2UV+), 
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where Vq = Lq^N2UU and V+ = (Lq - 2iLUo)-^N2UU. 
The amplitude-phase coupling factor S is obtained from 
the complex constant g and is given by 



(5 = 1^2 / 1^1 . 



(9) 



In this way, the factor 6 for STNO can be calculated 
numerically from the parameters of the LLGS equation. 

The noisy Hopf normal form given by Eq. ^ is de- 
rived when we add the noise term /(t) with {f{t)f{t')) — 
AD"f^5{t ~ t') to the SL Eq. (©. f{t) has the di- 
mension of frequency. The components in Eq. ([T]) are 
defined as W{t) = [D-f^ /vi)-^/'^W{t)e^^'^"+'^^-'^'^'^^^^\ 
i = y^D-f'^vit, p = F/-y/ D"f'^vi, and = u^/ D^f'^vi. 
Therefore, we can make the most of many well-known 
properties of Eq. ([1]!^°'— to examine the behavior of 
STNO. It is known, for example, that the spectrum 
linewidth Awfwhm far above the threshold {p ^ 0) is in- 
creased by a factor of (l-|-i5^)iiS In the context of STNO, 
when F 3> 0, the linewidth can be expressed as 



where = ^ det Lq . The right and left eigenvector of 
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FIG. 3: (Color online) (a) Dependence of 5 on the nonlinearity 
of the damping gi for various values of an external magnetic 
field H. An uniaxial anisotropy field is taken as Hk/4:TrMs = 
0.04. (b) Dependence of S on an external magnetic field H 
for various values of an uniaxial anisotropy field Hk- 



which corresponds to Eq. (11) in Rcf. 5. Here, fceT is 
the thermal energy. Awros is the linewidth at thermal 
equilibrium (luj — 0) given by Awrcs = 2rcq, where 



oq 



'T{loj = 0). Moreover, ii^osci is the magneti- 



zation oscillating energy and can be written as i?osci — 
2Un ^^%]f^ff^ U^oUPwVt,,, = h^^Pw when it is 
assumed that E'osci — near thermal equilibrium {en- 
ergy equipartition) . Here, Vfrcc is the volume of the free 
layer and Pw is the total power of W{t) given by Pw = 
a/D^^M{p + 2/F{p)} with F{p) = 0FeP'/4[l + erf(p/ 
2)]. From the expression of Eq. (fTO)l . it is found that 
the MR device in STNO itself is nothing but a resonator 
on the analogy of electrical circuits. The other one of 
well-known properties of Eq. ([1]) is that the amplitude- 
phase coupling factor distorts the power spectrum to non- 
Lorentzian lineshapes especially near the threshold (see, 
e.g., FIG. 5 of Ref. [HI). The degree of the hneshape 
distortion is determined by the magnitude of 6 and p, 
corresponding to the calculation in Ref. H. We com- 
ment on the validity of Eq. ([T]) for large-amplitude os- 
cillations. In Fig. [21 the theoretical fitting curves based 
on the model Eg. ([T]) are compared with the experimen- 
tal data of Ref. \1Q and give a good agreement with them 
up to / ~ 5.6 mA {p ~ 8.2) beyond the threshold cur- 
rent Ic = 4.8 mA {p — 0) estimated by the fitting. li 
Therefore, although the derivation of Eq. ^ is based on 
a perturbation expansion around the bifurcation point, 
it is considered to be valid for rather large-amplitude os- 
cillations with p 10. 

We briefly mention the oscillating frequency Wosci- 



From Eqs. ([T]) and ([7]), the oscillating frequency of a 
free layer magnetization far above threshold is written as 
Wosci — ujo^^S + rimXi. Although the calculation results 
for ImAi of Eq. ([5]) are not shown here, we have found 
that this quantity has a small value with ImAi ^ aQ 
for wide range of parameters of the LLGS equation. Ac- 
cordingly, Wosci is approximately given by tiJosci — loq — TS. 
Since F cx (I — Ic), while the frequency Wosci decreases as 
the current /(> Ic) increases when S > (red shift), Wosci 
increases when 6 < (blue shift) in accordance with the 
spin-wave models^i^i^i^. 

As illustrated above, the amplitude-phase coupling fac- 
tor S plays a key role to determine the behavior of an 
oscillator. Therefore, the features of STNO can be found 
out by the calculation of 6. 

Some calculation examples of S are shown in Fig. [3l 
It is considered the case where a free layer is an in- 
plane magnetic film with an in-plane external field ap- 
plied along the x direction, Jfcxt = Hx. It is assumed 
that TV = diag(0, 0, 1), ac = 0.02, and (0p, ^Jp) = (0, 0). 
In Fig. El^a), the dependence of S on the nonlinearity 
of the damping qi is shown. It is found that d mono- 
tonically decreases for qi and the variation of S is very 
large. This result suggests that a nonlinear damping sig- 
nificantly changes the LLC dynamicSfi^ In Fig.[3I^b), the 
dependence of 6 on an external magnetic field H for var- 
ious values of an uniaxial anisotropy field IIk{= ^K^^/ 
Ms) is shown. The nonlinearity of the damping is taken 
as qi — 3fi^ In the practical external field region, 6 is 
very sensitive to an uniaxial anisotropy field and varies 
largely. Therefore, when the dynamics of STNO is con- 
sidered, it is necessary to take the effect of an uniaxial 
anisotropy field into account seriously. This is the main 
result of the present paper. 

In summary, we have considered the dynamics of 
STNO by reducing the LLGS equation to a generic oscil- 
lator model and calculated explicitly the amplitude-phase 
coupling factor which is the key factor for the power 
spectrum. The amplitude-phase coupling factor S is very 
sensitive to magnetic fields, in-plane anisotropy, and the 
nonlinearity of damping. The large variation of S is the 
remarkable feature of STNO in comparison with conven- 
tional oscillators. The calculation way for S shown is ap- 
plicable for an arbitrary magnetization configuration and 
is useful for finding a stable STNO with small AtjpwHM 
(Eq. (flOl) ). which is preferable for applications. 
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The dimensionless power in Fig.[2][a) is given by P/RqI^ = 
a{p + 2/F{p)} with a ~ 2.3063 x 10"^ and p ~ 10.202(1 - 
/c). To obtain the linewidth in Fig.[2jb), we have used the 
parameters of sfTy^^l2-K — 11.24 MHz and S — 0.5, and 
have solved the eigenvalue problem of the Fokker-Planck 
equation corresponding to Eq. ([1} as done in Ref. [l3 or 
Ref.B 



